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§1.  Consider  the  supersonic  flow  over  a  "sphere -cone"  body 
(ABC  in  Fig.  1) ,  performing  plane  angular  oscillations  relative 
to  the  center  0  according  to  the  law 


a  =  Uq  cos  cjDt  (i-1) 

where  a  is  the  instantaneous  value  of  the  angle  of  attack.  In 
Fig.  1  ME  is  a  condensation  shock  and  QP  is  the  sonic  line. 

We  will  assume  that  the  conditions 


“O  «  ^ 


—  « 
Vi 


(1.2) 


are  satisfied,  where  Vp  is  the  free  stream  velocity  and  L  is 
a  characteristic  length  of  the  body.  Therefore,  the  perturbations 
caused  by  the  oscillating  body  will  be  small  so  that  the  problem 
can  be  solved  by  the  method  of  small  disturbances.  We  note  that 
the  second  condition  (1.2)  is  satisfied  to  a  high  degree  of  ac¬ 
curacy  for  supersonic  flight. 


Restricting  ourselves  to  the  linear  approximation  with  re¬ 


spect  to  the  frequency,  we  will  neglect  quantities  of  order 


a, 


0 


and  order  .  Then  the  parameters  of  the  gas:  the  velocity 

V,  the  pressure  p  and  the  density  p  can  be  represented  in 
the  following  form; 


V 


V, 


„  +  av 
0  a 


av. 

a 


p  =  Pq  +  ap^  +  ap^,  p  =  +  ap^  +  dp^  (1.3) 


The  parameters  with  index  0  describe  the  basic  field  aris¬ 
ing  from  the  stationary  flow  over  the  body.  Parameters  with 
indices  a  and  d  describe  the  perturbation  field  which  occurs 
in  phase  with  the  angle  of  attack  and  the  angular  velocity  respec¬ 
tively  and  do  not  depend  on  the  time  and  the  frequency. 

Let  us  substitute  the  expansions  (1.3)  into  the  system  of 
equations  of  gas  djmamics  written  in  a  moving  coordinate  system 
connected  to  the  body.  Then, for  determining  parameters  with  in¬ 
dex  a  we  obtain  the  system  of  linear  equations 
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grad(v«v  )  -  v„  x  rot  v  -  v  rot  v„  =  — ;r  grad  p„  - 

a  a  0  2°*^0 


—  grad  p 

Pq  ^ 


—  grad  (  —  -  y  — )  +  v  grad  (  — -  )  =  0 
T  0  ^  V,p_  p_y  a  ^ 


div(p„v  +  p  v„)  =  0 


and  for  determining  parameters  with  index  ct 


— >  ^  1  — >  ^ 

V  +  grad  v^(v.  -  v  .)  -  v„  x  rot  v. 
a  °  [  0  ^  a  ea''  0  a 


Pa  1 

"  "  ^ed>  ^  Pq  ■  j:  Pd 

Pq  "^0 


(1.4) 


(1.5) 


p  +  div(p.Vrt  +  PrvV.)  -  V  .  grad  p„  =  0 
0  *^0  a'  ea  °  *^0 


Pr>  P  P  -  P  • 


PfVx  — >  — >  Ppv 

Y  +  (v.  -  V  .)grad  —  =  0 

'  p  y  '  a  ea'*’  Y 

d  J  Pn 


where  Vg  is  the  moving  velocity  vector.  Here  the  parameters 
with  index  0  satisfy  the  system  of  equations 

o 


grad  -  Vq  X  rot  v^  = 


-  ^  Po 


(1.6) 


div  PqVq  =  0 


r\  ' 


which  describe  a  stationary  gas  flow. 
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We  will  suppose  that  the  spherical  and  conical  parts  of  the 
surface  are  joined  without  jumps  of  the  slope  generatrices  (with¬ 
out  corners)  .  Furthermore,  let  the  semiangle  of  the  cone,  0g, 
be  such  that  the  point  of  junction  B  is  always  found  in  the 
supersonic  region.  The  first  family  characteristic  BD,  start¬ 
ing  from  the  point  of  intersection  B,  divides  the  region  be¬ 
tween  the  condensation  shock  and  the  surface  of  the  body  into  two 
parts  for  the  stationary  flow  over  the  body  at  zero  angle  of  at¬ 
tack.  In  the  region  ABDM  the  flow  is  the  same  as  the  flow  over 

an  isolated  sphere  and  does  not  depend  upon  the  presence  of  the 
conical  part .  In  the  region  DBG  the  flow  is  determined  by  the 

motion  of  the  conical  part  of  the  surface  of  the  body  and  by  the 

propagation  of  the  initial  perturbations,  given  on  the  character¬ 
istic  BD,  and  also  by  the  stipulated  flow  over  the  spherical 
part  of  the  body. 

§2.  We  will  begin  with  the  determination  of  the  flow  in 
region  ABDM.  For  this  let  us  consider  the  flow  over  the  sphere 
performing  plane  apgular  oscillations  with  respect  to  the  center 
0  (see  Fig.  2)  according  to  the  law  (1.1).  The  motion  of  the 
sphere  can  be  decomposed  into  a  translational  motion  with  velocity 

v^  equal  to  the  absolute  velocity  of  the  center,  0]^,  of  the 

sphere  and  a  rotational  motion  around  the  center  with  velocity 
— 

v_  ,  determined  by  the  law  (1.1).  We  will  introduce,  in  the  ex- 
amination,  a  spherical  system  of  coordinates  R2,  ©2^  ^2  connec¬ 
ted  with  the  body  with  the  relative  velocity 


V  =  V  +  V 

e  e  e 


(2.1) 


and  a  partially^connected  system  of  coordinates  (spherical  co¬ 
ordinates  Rj^,  and  rectangular  coordinates  x,  y,  z) , 

translating  with  velocity  v^  (Fig.  2). 

Let  us  consider  the  perturbations  coinciding  in  phase, as  is 
defined, with  the  angle  of  attack.  In  the  linear  approximation 
with  respect  to  frequency  these  perturbations  correspond  to  the 
stationary  flovi?  over  a  body  at  an  angle  of  attack  a.  It  is 
evident  that  for  the  stationary  flow  over  the  sphere  the  gas  dy¬ 
namic  parameters  in  the  system  of  coordinates  Rj^,  9-^,  do  not 

depend  on  the  angle  of  attack  but  in  the  system  of  coordinates 
R2,  02-»  ^2’  linear  approximation  with  respect  to  a,  we 

have  the  relations 


Research  Department 
TR-24 

December  1962 


3 


f(a).fg+a(|£)  ,  v(a)-Vo+a(|j) 

a=0  a=i 


where  f  is  an  arbitrary  scalar  parameter.  Making  use  of  the 
relations  for  the  partial  derivatives  with  respect  to  time  in  the 
moving  system  of  coordinates,  we  obtain  for  the  scalar  parameter 
f 


f 

a 


V  . 

ea 


grad  f^ 


(2.2) 


and  for  the  vector  velocity  v 


V 


a-Cii) 


a=0 


=  (v  .  V)  v„ 
ea  0 


[k  X  Vq 


(2.3) 


or 


V 

a 


gradCVe^  Vq) 


V  .  X  rot  Vp. 
ea  0 


where  k  is  a  unit  vector  directed  along  the 

— 

tangular  system  of  coordinates,  and 


ea 


axis  of  the  rec- 
/d . 


Relations  (2.2),  (2.3)  permit  us  to  express  the  parameters 
with  indices  a  by  means  of  the  parameters  of  the  stationary 
flow  over  the  sphere . 

Let  us  now  consider  the  perturbations  coinciding  in  phase, 
as  defined,  with  the  angular  velocity.  The  rotational  motion  of 
the  sphere  around  its  center  does  not  cause  perturbations  in 

the  flow  of  an  ideal  gas.  Therefore,  the  perturbations  with 
indices  d  are  determined  by  the  translational  flow  over  the 
sphere  with  relative  velocity 


~^-k 

V 

e 


— > 

=  -  d^^j 


(2.4) 
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where  j  is  a  unit  vector  directed  along  the  y  axis  of  the 
rectangular  coordinate  system  and  £i  is  the  distance  between 
the  points  and  0,  measured  in  the  direction  of  the  axis 

©2=0  (see  Fig.  2) .  In  the  linear  approximation  with  respect 
to  frequency  the  translational  motion  of  the  body,  with  velocity 

(2.4)  in  the  partially  fixed  system  of  coordinates,  is  equivalent 
to  a  stationary  flow  at  a  fictitious  angle  of  attack 


a£^ 


(2.5) 


Taking  this  into  account,  we  obtain  for  the  arbitrary  scalar 
parameter  f 


f.  = 
a 


f 

a 


— ^ 

and  for  the  velocity  vector  v 


ea 


or 


— >  — > 

V.  -  V 

a 


ea 


V  . 


ea 


(2.6) 


(2.7) 


Relations  (2.6)  and  (2.7)  also  permit  the  parameters  with 
indices  d  to  be  expressed  by  means  of  parameters  of  the  sta¬ 
tionary  flow  over  the  sphere. 

It  is  easy  to  show  that  the  obtained  solutions  (2.3)  and 

(2.6),  (2.7)  satisfy  the  system  of  equations  (1.4),  (1.5)  and  the 
boundary  conditions  on  the  body 


aR 


=  0 


dR 


-  V 


eaR 


=  0 


for 


R  =  R 


(2.8) 


where  the  index  R  denotes  the  projection  of  a  vector  on  the  R 
axis  of  the  spherical  system.  In  fact,  because  the  system  of 
equations  (1.4)  and  the  first  relation  (2.8)  can  be  obtained  by 
variation,  with  respect  to  the  angle  of  attack  a,  of  the  exact 
system  of  differential  equations  of  gas  dynamics  and  the  boundary 
conditions  on  the  surface  of  the  body,  written  for  the  case  of 
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stationary  flow,  then  it  is  clear  that  the  solution  (2.2),  (2.3) 
satisfies  both  the  system  of  equations  (1.4)  and  correspondingly 
condition  (2.8).  Substituting  (2.2)  and  the  second  relation 
(2.3)  into  the  system  of  equations  (1.5)  and  using  relation  (2.1), 
we  obtain 


f  ,  1 

grad  v„ (v.  -  v  .) 

°  [  0  a  ea-^  J 


Vq  X  rot 


V  .)  X  rot  v„ 


grad  p  -  j-  grad  p 
Pq  "^0 


f  ^ 
div  (v. 

L  '^0  a 


— > 

-  V 


.) 

ea"^ 


p  .y 


OJ 


1=0 


(2.9) 


^0 

—  v„  grad  (  —  - 
y  0  ^  V^p^ 


y 

pQ 


P- 

PO'^ 


(v.  - 

^  a 


v.^)grad  -—  =  0 


ea' 


Comparing  (2.9)  with  (1.4)  and  the  second  relation  (2.8)  with 
the  first,  we  see  that  the  solution  (2.6),  (2,7)  satisfies  the 
system  of  equations  (2.5)  and  correspondingly  condition  (2.8). 

§3.  Let  us  now  dwell  on  the  calculation  of  the  flow  in  the 
region  DBG,  We  introduce  a  special  cylindrical  coordinate  sys¬ 
tem,  connected  with  the  body  (Fig.  1) .  The  axis  of  this  system 
coincides  with  the  axis  of  symmetry  of  the  body.  The  meridian 
plane  is  fixed  by  the  angle  ip  =  -  p.2,  and  the  position  of  a 
point  in  the  meridian  plane  is  fixed  by  polar  coordinates  R,  Q 
with  pole  at  the  point  of  intersection  of  the  generating  spheri¬ 
cal  and  conical  parts  of  the  body.  The  solution  of  the  system 
of  equations  (1.6),  (1.4)  and  (1.5)  must  satisfy  the  boundary 
conditions  on  the  surface  of  the  body 


''oa  ■  “  ’  ''ae  -  °  ’  ''da  '  ''eda  =  “  ®  -  »c 

where  the  index  G  denotes  the  projection  of  a  vector  on  the  9 
axis,  and  must  take  the  given  values  on  the  characteristic  surface 
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BD,  which  is  determined  as  a  result  of  the  calculation  of  the 
flow  over  the  spherical  part  of  the  body. 

Because  the  flow  in  the  region  DBG  is  supersonic,  then  for 
the  solution  of  the  nonlinear  system  of  equations  (1.6)  and  for 
the  solution  of  the  linear  system  (1.4)  and  (1.5)  we  may  use  the 
method  of  characteristics.  However,  this  method  for  obtaining 
numerical  results  requires  a  great  outlay  of  labor  and  the  appli¬ 
cation  of  high  speed  electronic  machines  in  each  concrete  case. 
Therefore,  we  cite  here  an  approximate  analytical  method  of  so¬ 
lution  of  the  problem.  A  solution  will  be  sought  in  the  forml 

f  =  (foo"*"  •  '  •) 

(3.2) 

+  ^^a^l  +  •  •  Osin  ^  +  d(v^^Q  +  rv.^3_  +  .  .  .)sin  ^ 


where  f  is  an  arbitrary  unknown  parameter,  excepting  the  projec¬ 
tion  of  the  vector  velocity  on  the  ^  axis.  The  coefficients  in 
the  expansion  (3.2)  depend  only  upon  0. 

Substituting  the  corresponding  expansions  into  the  system 
of  equations  (1.6),  (1.4)  and  (1.5),  written  beforehand  in  terms 

of  projections  on  the  axis  of  the  cylindrical  system  of  coordi¬ 
nates  R,  0,  ijj)  we  obtain  systems  of  ordinary  differential 
equations  with  respect  to  0.  Solutions  of  these  systems  must 
satisfy  the  boundary  condition  at  9=0^  and  the  initial  condi¬ 
tions  at  6  =  6*,  where  0*  is  the  angle  of  inclination  of  the 
tangent  to  the  characteristic  BD  at  the  point  B  (for  r  =  0)  . 
For  parameters  with  index  00  we  obtain  the  following  system  of 
nonlinear  equations 


This  form  of  dependence  upon  ijj  is  necessary  from  considera¬ 
tions  of  the  boundary  conditions  on  the  body  and  the  form  of 
the  initial  data  on  the  characteristic  BD  (see  the  following) . 
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'^OrO  '^oeo  " 


'  +  =  Pqo 

^000  '^OrO  Poo^oeo 


t  I 

Pqo  ^oeo  ^OrO 
Pqo  ^000 


0 


Pqo 

-  =  const. 


(3.3) 


and  the  supplementary  final  relation 

2  .  2 

^OrO  "^000  _y_  £oo  _ 

^  Pqo  ^ 

where  Vj^ax  the  velocity  of  the  flow  into  a  vacuum  and  the 

prime  denotes  differentiation  with  respect  to  0.  Let  us  substi¬ 
tute  the  expressions  for  Pqq/PoO’  PoO^PqO  POO'^POO 

equations  (3.3)  and  (3.4)  into  the  result  of  differentiating 
(3.4).  Taking  into  account  the  first  equation  (3.3),  we  obtain 


max 


/"j  /,  \ 


''orO  ''OrO  '  ''ma^^''0r0  +  ''0r0>  ‘  ° 

Equating  to  zero  the  first  factor,  with  the  aid  of  (3.3) 
and  (3.4),  we  obtain  a  Prandtl-Meyer  solution.  This  solution  for 
parameters  with  index  00  would  be  required  by  us  in  the  case  of 
a  junction  of  the  spherical  and  conical  sections  with  a  jump  in 
slope  at  the  point  B.  In  the  case  under  consideration  it  is 
necessary  to  equate  to  0  the  second  factor  (3.5).  Thereby  we 
obtain  a  solution  satisfying  the  initial  and  boundary  conditions: 


V, 


V 


* 
OrO 


OrO 


cos  cp. 


cos  9 


'^000 


'OrO 


cos  9 


srn  9 


Poo 


(9=0-  ©c) 


(3.6) 
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where  the  asterisk  superscript  denotes  values  of  the  parameters 
at  0=0*.  The  solution  (3.6)  describes  a  uniform  translational 
motion.  For  the  remaining  coefficients  in  the  expansion  (3.2)  we 
obtain  a  system  of  linear  differential  equations  which  are  easy 
to  integrate  in  finite  form.  We  will  cite  the  results  of  the 
integrations  for  the  coefficients  in  the  expansion  of  the  pressure 
p,  necessary  for  the  calculation  of  the  aerodynamic  loads  acting 
on  the  body 

*  * 

PqI  *  Pal  * 

PqI - i  PqO  =  PqO^  Pal  =  - ^  PdO  =  PdO 

cos  cp  cos  9 

(3.7) 

9 


dl 


cos  9 


■^fp 


dl  ^PQO'^OrO 


*  tan  9*] cos  9 


-  2p, 


OrO 


* 


00  * 
cos  9 


sin 


For  the  completion  of  the  solution  of  the  problem  in  region 
DBG,  it  is  necessary  to  express  the  initial  data  on  0=0* 
with  the  aid  of  the  solution  obtained  in  §2  for  the  flow  in  re- 
gi.on  ABDM. 

§4.  Let  us  begin  with  the  derivation  of  the  equation  of  the 
actual  perturbed  characteristic  surface,  passing  through  the  line 
of  intersection  of  the  spherical  and  conical  parts  of  the  surface 
of  the  body.  We  will  seek  this  in  the  form: 

"fC 

9  =  (_G  +r0j^+...)  +  +  r0^j^  +  .  .  .)  cos  f  +6(0^^  +  r0^^  +  .  .  . )  cos  ^(4.1) 


The  form  of  the  characteristic  surface  is  determined  by  the 
condition 


N-n(v-v)=a, 

e 


y  I  grad  F|2 


n 


grad  F 
y  I  grad  F  |  2 


(4.2) 


where  n  is  the  unit  normal  vector,  N  is  the  relative  velocity 
of  motion  of  a  point  of  the  surface  in  the  direction  of  the  normal. 
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F ir , 6 , f , t)  =0  is  the  equation  of  the  characteristic  surface 
and  a  is  the  local  velocity  of  sound. 


Substituting  (4.1)  into  the  second  and  third  relations 
(4.2),  we  obtain 


N  =  -  d(re^Q  +  ...)cos  ^ 

=  -  (r0^  +  . .  .)  -  a(r^]^  +  . .  .)  cos  -ip  -  a(r0^^  +  .  .  .)cos  f  (4.3) 


’^S  ^-ip  ^  R —  ^  —  +  .  .  .)  srn  Ip 


where  Rq  =  BF  (Fig.  1).  Let  us  substitute  the  expansions  (3.2) 
and  (4.3)  into  the  first  relation  (4.2).  Expanding  the  functions 
of  0  in  a  series  in  terms  of  0-6*  by  (4.1),  we  obtain  ex¬ 
pressions  of  the  coefficients  in  the  expansion  of  the  equation 
of  the  characteristic  surface  (4.1)  in  terms  of  the  values  of 
the  parameters  of  the  gas  at  6  =  6*.  We  cite  the  expressions 
for  two  of  the  coefficients  which  are  necessary  to  us  in  what 
follows 


*  •* 
^a60  ^aO 
aO  “  “  '  *  ,  ’ 

’^oeo  ^00 


^deo 


+  a 


do 


dO 


V, 


060 


•  * 
^00 


(4.4) 


^  r\>  3.re  the  coefficients 

00  aO  aO 

velocity  of  sound  (3.2). 


in  the  expansion  of  the  local 


Writing  down  the  condition  of  continuity  of  parameters  on 
the  actual  perturbed  characteristic  and  expanding  the  functions 
of  6  in  a  series  with  respect  to  6  -  6*,  we  obtain  that  at 
6=6*  it  is  necessary  to  satisfy  the  following  relations: 
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0 


00 


=  0, 


[f  * 

[^01 


[  I*] 


Ko]  +  ‘’ao'Cj  -  0 


\f  1  +  U'*]  +  o  if'*!  +  B  o  ff"*!  +  if'*!  =  0 

i^alj  ^  aOj  ^  ^aoi^OlJ  ^^l^aoi^OOj  ^  ®al[^00j  ° 


(4.5) 


ff  *1  .  0  ff'*l  =  0 

i  aOj  ^  ®doi^00j  ^ 


r  *1  f'*l  f'*l  f^'*l 

[^dlj  ®l[^dOj  ^dOL^OlJ  ■*■  ®l®d0[^00j  ■*■  ®dl[^00j 


=  0 


where  the  bracket  denotes  a  jump  in  the  quantities  contained  in 
it  from  the  left  to  the  right  for  the  approximation  at  0=0*. 
The  continuity  of  the  parameters  of  the  gas  on  the  actual  per¬ 
turbed  characteristic  surface  and  the  distinct  analytical  char¬ 
acter  of  the  solution  (jumps  in  derivatives)  on  both  sides  of  it 
lead  to  the  fact  that  for  a  united  solution,  describing  the  flow 
on  the  spherical  and  conical  parts  of  the  body  it  is  necessary 
to  take  into  account  discontinuities  of  parameters  along  the  line 
0  =  0*,  as  indicated  in  (4.5). 

In  the  conclusion  of  this  paragraph  let  us  make  one  remark. 
If  the  system  of  equations  (1.6),  (1.4)  and  (1.5)  is  projected 
on  the  axis  of  the  cylindrical  system  of  coordinates  and  if  one 
eliminates  from  the  parameters  of  the  perturbations  (with  indices 
a  and  d)  the  dependence  upon  ip,  then  it  is  easy  to  see,  that 
the  coefficients  for  the  corresponding  derivatives  in  all  three 
systems  coincide,  and  consequently,  the  characteristics  of  the 
linear  systems  of  equations  for  the  perturbations  coincide  in 
the  physical  plane  with  the  characteristics  of  the  nonlinear 
systems  of  equations  describing  the  stationary  flow  over  a  body 
at  zero  angle  of  attack.  This  circumstance  is  convenient  for 
the  integration  of  the  linear  system  of  equations  (1.4)  and  (1.5) 
by  the  method  of  characteristics.  However,  from  the  preceding 
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consideration£5 j  it  is  clear  that  for  calculations  with  the  aid  of 
relations  of  the  type  (4.5)  it  is  necessary  to  take  into  account 
discontinuities  of  the  perturbation  parameters  on  the  character¬ 
istics  separating  regions  with  different  analytical  character  of 
the  solutions . 

§5.  With  the  results  of  §2,  the  coefficient  values  in  the 
expansion  (3.2)  are  determined  on  the  line  9  =  d* ,  by  using  the 
approximation  on  the  side  of  the  spherical  part  of  the  body. 

These  coefficient  values  will  be  denoted  by  the  subscript  s, 

^In  the  vicinity  of  the  surface  6  =  6"*^ ,  9^  »  cl,  and 
TT  “  ©2  »  so  that  |j.2  l-‘-l  4-  0(a)  .  Taking  this  into  account, 

and  remembering  the  expressions 

**  ic-k 

v.„=0,  v._=-R  cos  p  ,  V  ,  =  R  cos  9  sin  p  (5.1) 

eaR  ^  ea9  e.af  ^  ' 

from  (2.3)  and  (2.2)  we  find 


'V,  r\_,  U  J-Q 

f  ^R2 j  ^0 ^^2 ^  ^2^  *  ^  cos  ^2^ 

d02 

~  otVog 

V  (R2,  ^2’  ^2’  ^  ^  ^2 

^  sin  ©2 


(5.2) 


One  can  easily  find  the  desired  values  of  the  coefficients 
for  the  expansion  (3.2)  for  0=0*  on  the  side  of  the  spherical 
part  with  equations  (5.2)  and  (2.6)  and  using  the  relations  between 
components  of  the  velocity  vector  in  cylindrical  and  spherical  co¬ 
ordinate  systems 


v^  =  -  Vg^  cos(0  +  9^)  +  Vg  sin(e  +  9^) 

(5.3) 

Vg  =  sin(0  +  ©2)  +  Vg  cos(0  +  9^)^7^  =  -  v^ 
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and  the  relations 


C^)  -  'I'  . 

r=0 


C— )  „  -  0 

r=0 


(5.4) 


Ca?^ 


COS  cp 

R 


where  R  is  the  radius  of  the  sphere. 


Assuming  R  =  1  and  using  the  relation  (69)  ^ 

o  ^  ir=U 

for  coefficients  in  the  pressure  expansion,  we  have 


=  iir  -  Oc., 


^PqO^s  ^P0^r=0 


(Poi)  ^  -  sin  <P  (;^)  +  cos  T 

2  r-0  d02  r=0 


(Pqi)  -  00s  <p 


.  * 


.  3ln 

2'^r=0  be/ 


2  r=0 


<PaO>s  -  -  /  J 

2  r=0 


(5.5) 


-  srn  9 


■*  / 

-  cos  9  ( 

2^'^2  r=0 


2  r=0 


(PdO^s  - 


f 

V^'^PqO-'c  ' 


f  _  1/  N* 
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Equations  (4-5)  determine  the  relation  between  values  of  the 
coefficients  of  expansion  in  r  on  the  line  0=9*  by  approach¬ 
ing  the  line  from  both  sides.  For  the  pressure,  in  particular, 
we  have 


P00*<P00>s’  Pfll"  <P01>s'  PaO’  (PaO^s’  Pal  "  ':Pal>  s ®a0  |Poi] 


*  *  *  * 

PdO  *  <PdO>s’  Pdl  ■  (Pdl^s  +  ®d0lP01J 


r  1*1 

The  magnitude  of  the  jump  in  [PoiJ  these  formulas, 

according  to  (3 . 7)  is  determined  by  the  relation 


(5.6) 


[1*1  1*  *  * 

Poll  ■  <Poi>s  +  »  (Poi^s 


(5.7) 


and  the  parameters  G  „  and  0.^  by  the  relations  (4.4) 


^00 


=  -  cos  Cp  (Vq~)^^q 


=  -C^)  ' 

^^2  r=0 


^^2  r=0 


(5.8) 


^dO 


_i  * 


-  ^  ^aOO  +  ^ 


§6 .  Let  us  consider  how  to  determine  the  aerodynamic  moment 
M2,  acting  on  the  oscillating  body  in  supersonic  flow.  The 
moment  is  made  up  of  two  parts.  The  first  part  (Mz)s  is  crea¬ 
ted  by  the  pressure  force,  acting  on  the  spherical  part  of  the 
body  surface,  the  second  is  the  pressure  force  acting  on  the 
conical  part  of  the  body  surface.  We  introduce  the  pressure  co¬ 
efficient  by  the  formula 
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(6.1) 


M 


C 

mz 


5PlV^LTr(L  tan  9^  ) 


where  is  the  density  in  the  undisturbed  flow^  and  L  is  a 

characteristic  length  of  the  body  (Fig.  1) . 

Then  the  derivatives  of  the  moment  coefficient  with  respect 
to  a  and  P,  where  P  =  otL/V]^,  can  be  represented  for  the 
spherical  part  in  the  form 


^^mza^s 


=  4 


1  + 


X 


(7-1)M^J 


-i-x 


tt/ 2-6 


cos  9 


(Po)r  _Qjsin  9  cos  9  d0 


2 

cos  6, 


(6.2) 


^^mzp^  s 


-a- 


X 


-0 

L 


2  o 

cos  9 


-)ic  ) 

y  ^  mza' s 


and  for  the  conical  part  in  the  form 


(G  )  =  2 

mza'  c 


1  + 


(Y-1)M[ 


a, 


0 


sin  9  cos  9 


a. 


.  2 

sxn  9 


rA  ^  ^ 

2V.L"  L  ^  ®c) 


(6.3) 
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For  convenience  of  comparison  with  the  aerodynamic  charac¬ 
teristics  of  conical  bodies  the  geometrical  parameters  of  the  cir¬ 
cumscribed  cone  are  introduced  in  formulas  (6.2)  and  (6.3)  (see 
Fig.  1)  . 

The  method  being  suggested  was  used  in  the  calculation  for 
=  4.0.  The  necessary  data  on  the  parameters  of  the  stationary 
stream  line  of  the  sphere  are  taken  from  paper  [IJ.  In  Figs.  3 
and  4  are  presented  the  functions 


C  =  (C  )  (C  ) 
mza  mza^  s  '  mza' 


c^ 


C  ft  =  (C  o)  + 
mzp  mzP  s 


c 


for 


BC  =  R 


(P  =  dL/V^) 


of  the  cone  half  angle  0^,  and  of  the  center  of  oscillations 
i/L.  A  positive  sign  on  0^7.^  corresponds  to  statistical  stabil- 
ity^  and  a  positive  sign  on  Cmzp  to  damped  oscillations.  The 
graph  in  Fig.  5,  where  the  value  of  Cj^zp  x  10^  is  given,  shows 
the  contribution  of  the  spherical  part  of  the  body  to  the  aero¬ 
dynamic  moment.  This  contribution  is  small  for  large  9^,  but 
rapidly  gtows  with  decreased  angle  of  attack,  if  only  the  center 
of  oscillation  is  not  too  closely  located  to  the  center  of  the 
sphere . 

§7.  Now  we  will  consider  the  flow  over  the  sphere,  per¬ 
forming  plane  angular  oscillations  with  center  0  (see  Fig.  2) . 
Tile  oscillations  about  a  certain  fixed  angle-of-attack  a*  are 
given  by  the  law 


a  =  a  +  AQq  cos  cot  (7.1) 

In  accord  with  Part  2  the  perturbations  caused  in  the  flow 
by  the  angular  oscillation  of  the  sphere  with  center  0  are 
equivalent  to  translation  with  absolute  velocity  by  the  sphere 
with  center  In  the  case  of  slow  oscillations,  if  one  is 

restricted  to  a  linear  approximation  with  respect  to  frequency, 
translations  of  the  body  cause  the  same  perturbations  in  the  gas 
flow  as  an  imaginary  change  in  both  angle-of-attack  and  velocity 
of  the  undisturbed  flow: 
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^  cos  a  ^ 

Aa^  =  -  - - -  a,  AV,  .  =  i,  sin  a  d 

$  1<1>  1 


(7.2) 


where  is  the  distance  between  points  0  and  O-^,  positive 
at  the  rear  position  of  the  center  of  oscillations . 

The  moment  coefficient,  relative  to  the  point  0,  aq^ting  on 
the  stationary  streamline  of  the  sphere  under  angle-of-attack  a, 
is  determined  by  the  relation 


s 


(7.3) 


where  the  radius  of  the  sphere  is  considered  as  the  characteris¬ 
tic  dimension,  Cjj  is  the  drag  coefficient  of  the  sphere.  Cal¬ 
culating  the  change  in  moment  caused  by  a  fictitious  change  of 
angle-of-attack  and  flow  velocity  (7.2)  we  obtain 


ac  1 


■'mz  V , 


c- ^  Cr:) 


or  taking  into  account  (7.3), 


1  * 

C  =  -  —  cos  a  C 
mza  R  X 

s 


dR 


(7.4) 


£  I  ? 

<=mz|3  -  (r-)  <=x  +  “l  ^ 
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In  conclusion  we  note  that  the  preceding  derivatives  of  the 
moment  coefficient  are  associated  with  the  generally  used  coeffi- 

a  B 

cients  of  rotational  derivatives  C  ,  C  C  in  the  follow- 

mz  mz  mz 

ing  manner  (see  [2j); 


C 


mza 


=  -  G 


a 

mz 


^mzB 


=  C  ^  + 
mz 


mz 


(B  =  a  L/V) 


(7.5) 
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